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1 Abstract 



^**^ ' The Refined Kolmogorov Similarity Hypothesis is a valuable tool for the description of intermittency in isotropic 



conditions. For flows in presence of a substantial mean shear, the nature of intermittency changes since the process of 
energy transfer is affected by the turbulent kinetic energy production associated with the Reynolds stresses. In these 
^ , conditions a new form of refined similarity law has been found able to describe the increased level of intermittency 
~j ' which characterizes shear dominated flows. Ideally a length scale associated with the mean shear separates the two 
ranges, i.e. the classical Kolmogorov-like inertial range, below, and the shear dominated range, above. However, 
the data analyzed in previous papers correspond to conditions where the two scaling regimes can only be observed 
, individually. 

' In the present letter we give evidence of the coexistence of the two regimes and support the conjecture that 

', the statistical properties of the dissipation field are practically insensible to the mean shear. This allows for a 
theoretical prediction of the scaling exponents of structure functions in the shear dominated range based on the 
known intermittency corrections for isotropic flows. The prediction is found to closely match the available numerical 
' and experimental data. 

o . 

C5 2 Introduction 

^ At large Reynolds number {Re) turbulent flows are characterized by strong non Gaussian intermittent fluctuations. 
^ For homogeneous isotropic turbulence, a quantitative measure of intermittency can be given by using the structure 
functions < 5V'P{r) > where 



■ SV{f) = [u{x + r) - u{x)] (1) 



r 



Then, the generalized dimcnsionless flatness 



F (r) < > (2) 



exhibits intermittency, in the sense that Fpir) — ^ oo for r — > and Re —t oo. For rj <ti r <^ Lq, where rj is the 
Kolmogorov dissipation length and Lq is the integral scale of turbulence, structure functions show scaling behavior, 
i.e. < 6VP{r) >cx r^'-P\ where C{p) are anomalous scaling exponents (C(p) 7^ p/lCil)), and C(3) = 1 due to the 
Kolmogorov four-fifth equation Q. It is a remarkable result, obtained in the last ten years, that C(p) are observed to 
be universal, i.e. independent of the Reynolds number and of the forcing mechanism, for homogeneous and isotropic 
turbulence 1^. 

Much less information is available for non isotropic turbulence. Recently, a number of experimental and 
numerical investigations [Q in shear flow turbulence have shown that intermittency increases when the shear strongly 
affects the energy cascade. In the language of scaling exponents, an increase of intermittency means that Cip) are 
different from those observed in homogeneous and isotropic turbulence. Based on DNS of turbulent channel ffow, it 
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was recently proposed Q that the increase of intermittency is due to the breaking of the Kolmogorov Refined Similarity 
Hypothesis (RKSH), which for homogeneous and isotropic turbulence reads |^ 

< SVPir) >oc< > r^/^ (3) 

where 

1 r 

eioc{x)(fx (4) 



1 

W) 



B(r) 



and B(r) is a volume of characteristic size r while eioc{x) is the instantaneous local rate of energy dissipation. Equation 
(^), in its Extended Self Similarity (ESS) formulation 

has been successfully checked for a wide range of Reynolds numbers in different homogeneous and isotropic turbulent 
flows 0. 

For shear flow turbulence, it has been suggested that equation (||) should be replaced by 

< SVPir) >oc " < SV^ >P/2 (6) 



for r ^ Ls where Lg — y^e/S^, S being the mean shear in the system. Following these considerations, the shear 
scale ideally separates, with regard to the scaling properties of structure functions, the range of scales where turbulent 
kinetic energy production prevails (r 3> is) from the range of scales characterized by purely inertial energy transfer 
(r ^ Ls). Equation is able to explain the observed increase in intermittency, by assuming that the scaling 
properties of < ej? > are not changed by the mean shear S. 

It is the aim of this letter to present and discuss numerical and experimental evidence that the two scaling regimes, 
predicted by eq. (|^) and (j^), are indeed observed simultaneously in the range below and above Lg respectively. 
Moreover, the results confirm the conjecture that the scaling properties of the energy dissipation are not modified 
appreciably with respect to homogeneous and isotropic turbulence. 



3 Data sets and analysis 

We discuss two sets of data, one obtained by a long and highly resolved DNS of homogeneous shear flow turbulence 
1^ , the other by hot wire measurements in the logarithmic region of a turbulent boundary layer ||] . 

Concerning the homogeneous shear flow, we have considered a turbulent flow with an imposed mean velocity 
gradient S free from boundaries. The Navier-Stokes equations, written in terms of velocity fluctuations, are solved by 
using an efficient pseudo-spectral method with a third order Runge-Kutta scheme for time advancement following the 
transformation of variable proposed by Rogallo [p^ . As shown by Pumir |TT| and recently confirmed by the present 
authors the flow reaches a statistical steady state characterized by large fluctuations of the turbulent kinetic 
energy. The growth of turbulent kinetic energy is associated to large values of the Reynolds stresses, produced by a 
well defined system of streamwise vortices via a lift-up mechanism [|l2|. In this flow, because of shear scale fluctuations 
due to the mentioned behavior of both turbulent kinetic energy and Reynolds stresses, the crossover between the two 
scaling ranges is not sharply defined. In fact, we observe an overlapping of the two scaling regimes, and the resulting 
scaling shows an effective slope. In order to reduce as much as possible the fluctuations of the shear scale, a conditional 
sampling is introduced by considering only flow configurations where the production term exceeds a given threshold. 
Among these configurations, only those corresponding to a large value of turbulent kinetic energy [8 > aSrms) are 
retained to reduce the mean value of the ratio Ls/rj cx fl^^^ (fl being the mean enstrophy). 

Concerning experiments, we analyze the velocity data on a flat plate boundary layer measured in a wind tunnel 
(test section length of 150cm) operated at 11.9to/s. The boundary layer thickness is ~ 25mm and the Reynolds 
number based on the momentum thickness is about 2200. The boundary layer has the expected logarithmic region 
with the usual log-law constants p3| . Hot wire measurements were performed at several distances from the plate, 
using a constant temperature anemometer. The data acquisition was long enough to achieve convergence of the sixth 
order structure function. 



2 



4 Double scaling regime 



We begin by analyzing the DNS of the homogeneous shear flow. We have strong evidence that for r > i^, eq. (Q) fails 
and the new form of RKSH is established, as reported in and further discussed in the following. Furthermore, 

the statistical properties of energy dissipation < > are not distinguishable from those observed in homogeneous and 
isotropic turbulence. The last statement can be directly checked by looking at figure [| where we plot < > versus 

< > both for homogeneous shear flow and isotropic turbulence while, in the inset of the same figure, we plot < > 
versus r/ry for both cases. At variance with DNS, a direct measurement of < > is not available for the experimental 
data and we are not fully confident in the one dimensional surrogate of ejoc as a direct measure of the local rate of 
energy dissipation, being the fiow strongly anisotropic. Nevertheless by using the one dimensional surrogate, we can 
practically reproduce the results shown in figure |l|. 

At any rate, to be cautious, we may avoid the explicit use of the energy dissipation by plotting structure functions 
in the form suggested by Ruiz-Chavarria et. al. p^ . Specifically, here, we introduce indicators based on < SV^lr) > 
to detect the two scaling regions and to compare our findings with the predictions made in eq. and (^). We remark 
that, both for numerical and experimental data, the Reynolds number is not large enough to observe the scaling of 

< SVP{r) > and < ef. > with respect to separation. Thus we employ the ESS as a valuable tool to estimate the scaling 
exponents. In particular this implies that the exponents r(g) are defined by the relation < >oc< 5V^ >'^^'^\ 

Following eq. (||) and (^) and the above discussion, we compute both for the DNS and the experimental data the 
quantity =< SV^ > / < 5V'^ >p^^ and pp =< SV^ > / < SV^ >p/^ which are expected to satisfy the relations: 

ap oc <^ (7) 

[ < ^-^^3 >r(p/3)-r(2/3)p/2 ^ ^ 

and 

r < ^-1/3 >r(p/2)~r(3/2)p/3 ^ ^ 

Ppoc < (8) 
[ < SV^ >^(p/3) r < L,. 

Equations (|^) and (|^) allows us to compare the ESS exponents of ap and pp against the exponents predicted by 
equations (||) and (^. Let us remark that eq. and (||) are also based on the assumption that T{q) are the same 
both for shear dominated fiows and homogeneous and isotropic turbulence, as supported by the DNS data for < > 
shown in figure |l]. 

In figure |^ we plot log erg against log < SV'^ > for the data of the homogeneous shear fiow and the turbulent 
boundary layer. The fits for r 3> are in close agreement with the value of r(3) = —0.59 expected from homogeneous 
and isotropic turbulence. In the inset of figure |^ we show the local slope d[logcr6]/'^[log < SV'^ >] computed from the 
homogeneous shear flow dataset. The two dashed lines indicate the numbers r(2) — 3r(2/3) and r(3), i.e. the expected 
scaling exponents for r <^ Ls and r ^ Lg respectively. 

Figure ^ shows the main result of this letter, i.e. the clear evidence of the two scaling regions below and above 
Lg. Specifically, the experimental and numerical evidence of the coexistence of two different intermittent regions in 
shear flow turbulence. As a further check of the theory, in figure ^ we show log pg versus log < 6V^ > both for the 
homogeneous shear fiow and the turbulent boundary layer while in the inset we show the local slope (i[log p6]/f^[log < 
5V^ >] . Also for the variable pg we can claim a very good agreement between the observed experimental and numerical 
results against theoretical predictions. Our results concerning the double scaling of structure functions are consistent 
and complementary with the generalized structure function < {5V^ + arSV'^y/^ > proposed by Toschi ct. al. [p^ . 

Finally we remark that, by denoting C(p) the anomalous exponents in shear flow turbulence and (^(p) the anomalous 
exponent in homogeneous isotropic turbulence, the ESS estimate becomes 



2 M 2 



The first term comes from the fact that in shear turbulence, we may express 

< SV"" >0C< £3/2 X jy3 >3C(2)/2 (^0) 
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which implies ^(2) = 2/3[l — t(3/2)]. Equation provides a theoretical estimation for the scaling exponents of 
structure functions in shear dominated flows by using the intermittency corrections of isotropic turbulence. The 
values given by eq. (^) are compared against their direct measure in the DNS of the homogeneous shear flow in table 

i 

5 Final remarks 

We have analyzed scaling properties of velocity fluctuations in shear dominated flows considering a DNS of a homoge- 
neous shear flow and hot wire measurements in the logarithmic region of a turbulent boundary layer. When dealing 
with DNS data a conditional statistical analysis has allowed us to consider only shear dominated samples and to reduce 
the fluctuations of the shear scale resulting in a well defined scaling behavior below and above Lg respectively. To by 
pass the use of the one dimensional surrogate of energy dissipation, in the experiments, we have evaluated the refined 
similarity laws only in terms of structure functions. 

The analyzed data show clearly the coexistence of a double scaling behavior of structure functions across the shear 
scale. Actually the scaling exponents computed below and above the shear scale are in a very good agreement with the 
theoretical values provided by the classical RKSH |^ , , and the new refined similarity law respectively. Since the 
statistical properties of the energy dissipation are weakly affected by the shear, intermittency corrections of isotropic 
fiows can be successfully employed in the new form of similarity law, see table |^, to predict the scaling exponents of 
structure functions in the shear dominated range. 
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Table 1: Scaling exponents of structure functions (DNS) above and below the shear scale Lg- Data are compared 
with those of homogeneous and isotropic turbulence and with the prediction of eq. (H) . 



p 


1 


2 


3 


4 


5 


6 


r < Ls 


0.36 


0.69 


1.00 


1.28 


1.54 


1.78 


r>Ls 


0.38 


0.72 


1.00 


1.23 


1.42 


1.56 


hom.iso 


0.36 


0.69 


1.00 


1.28 


1.54 


1.78 


eq. (PI) 


0.39 


0.73 


1.00 


1.23 


1.42 


1.58 
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Figure 1: < > vs. < > in the homogeneous shear flow (circles) and in homogeneous isotropic turbulence 
(diamonds). In the inset < > vs. r/rj for the two cases. 




Figure 2: logCTg vs. log < SV^ > in the homogeneous shear flow (circles) and in the turbulent boundary layer at 
y"*" = 115 (triangles). DNS and experimental data are fitted at scales r > Lg by power laws with a slope s = —0.58 and 
s = —0.59 respectively. In the inset, the local slope (i[log(T6]/d[log < 5V^ >] vs. r/rj in the homogeneous shear flow 
obtained by considering the conditional sampling with a = 1.3 (solid line). The dotted lines corresponds to the two 
scalings given by eq. (|^) at scales r < Ls (—0.3) and r > Lg (—0.59) using the values of T(g) for isotropic turbulence. 



6 



4 6 



Figure 3: logpg vs. log < SV^ > in the homogeneous shear flow (circles) and in the turbulent boundary layer 
at y+ = 115 (triangles). Both DNS and experimental data are fitted at scales r < Lg hy a, power law with a slope 
s = —0.22. In the inset, the local slope d[\og pe] / d[\og < SV^ >] vs. r/rj for the homogeneous shear flow obtained 
by considering the conditional sampling with a = 1.3 (solid line). The dotted lines corresponds to the two scaling 
given by eq. (||) at scales r < Lg (—0.22) and r > Lg (—0.41) using the values of r(g) for isotropic turbulence. 
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